Modification of electromagnetic structure functions for the yZ-box diagram 
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The yZ-box diagram for parity violating elastic e-p scattering has recently undergone a thorough analysis by 
several research groups. Though all now agree on the analytic form of the diagram, the numerical results differ 
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due to the treatment of the structure functions, F l 2 3 (x , Q ) . Currently, F l2 3 { x , Q ) a t low Q and W must be 
approximated through the modification of existing fits to electromagnetic structure function data. We motivate 
and describe the modification used to obtain Fj (x, Q ) in our previous work. We also describe an alternative 
modification and compare the result to our original calculation. Finally, we present a new modification proce- 
dure to acquire Fj {x,Q ) in the resonance region and calculate the axial contribution to the yZ-box diagram. 
Details of these modifications will illuminate where discrepancies between the groups arise and where future 
improvements can be made. 



I. INTRODUCTION 

Parity violating e-p scattering experiments performed at 
momentum transfers away from the Z-pole are used to test 
the Standard Model prediction of the running of sin 9w- 
The Q we ak experiment at Jefferson Lab [1] aims to perform 
a 0.3% measurement of sin 9w at a momentum transfer of 
Q 2 = 0.026 GeV~. To obtain this desired precision, all radia- 
tive corrections must be known to an even higher precision. 
Up to one loop order, the weak charge of the proton at zero 
momentum transfer is given by H 

Q^ = (1+Ap + A e )(l -4sin 2 0w(O)+A;) 

+ CW + n zz + ReD yZ . (1) 

Here, A e and A' e are electron vertex corrections, Ap is the 
W and Z mass renormalization, and 1—4 sin 9w (0) is the one 
loop value of the weak mixing angle evaluated at Q 2 = 0. The 
WW and ZZ box diagrams, Oww and Dzz, are dominated by 
large momentum exchange and can be calculated using pertur- 
bative QCD. A different technique is required to calculate the 
yZ-box diagram due to low Q 2 contributions. Gorchtein and 
Horowitz J3I] used a dispersion relation to evaluate the yZ-box 
diagram at zero momentum transfer and obtained a result that 
was larger than expected B- Sibirtsev et al. |01 used the same 
technique and found an analytic result that was greater by a 
factor of 2. This discrepancy inspired a third calculation |J5[] 
that agreed with the Sibirtsev et al. result. After reevaluating 
their work, Gorchtein et al. fla] confirmed the factor of 2. All 
three groups now agree on the analytic form of the yZ-box. 
The imaginary vector portion is 
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where 
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A = 



(2ME) Z - 2ME(W Z -M 2 + Q 1 ) -M l Q 
Q 2 {W 2 -M 2 + Q 2 ) 



(2) 



(3) 



In the above equations M is the mass of the proton, E is 
the lab energy of the incoming electron, s = M~ + 2ME, 



W% = (M + m K ) , m K is the mass of the pion, and Qf„ ax = 
(s -M 2 )(s - W 2 )/s. The fine structure constant a em (Q 2 = 0) 
is used because the integral receives most of its support from 
low Q 2 . The dispersion relation that relates Im OY, z to Re nL 

is 



ReD v rZ (E) = 



IE 



dE' 



E' 2 -E 



jlmD v rZ (E') 



(4) 



where v„ = {W 2 -M 2 )/2M. 

The Q wea k experiment ran at an incoming electron energy 
of E = 1 . 165 GeV. TableQ]shows the numerical Re OY, z results 
obtained by each group at this energy and includes a soon to 
be published value 121]. The differences occur because of the 
models used for the Fj structure functions. Currently, there 
are no data for these structure functions at low Q 2 and W 2 
and each group performed calculations using their own mod- 
ifications to electromagnetic structure functions. The PVDIS 
experiment [SO at Jefferson Lab has several data points for the 
deuteron's FJ 2 3 m the resonance region. These data will be 
insufficient to produce a model-independent fit, but provide a 
first step in testing the validity of the modifications |9||. 



TABLE I: ReD^ z x 10 3 evaluated at E = 1.165 GeV. 



Sibirtsev et al. [4J 4.7+q \ 

Rislow and Carlson [5J 5.7 ± 0.9 

Gorchtein et al. [6J 5.4 ± 2.0 

Hall et al. \T\ 5.60 ±0.36 



The axial contribution to the yZ-box has also recently un- 
dergone analysis. The axial contribution to ImD r z is 



rZV ; {2ME) 2 J M i 



dW l 



QlmX d0 2 a (o 2 )Sv(Q 2 ) BF f^Q 2 ) 
dQa em (Q) - l+Q2/M 2 



where 



B = 



2ME 



W 2 -M 2 + Q 2 2 



(5) 



(6) 



TABLE II: ReD^ z x 10 3 evaluated at E = 1 . 165 GeV. 



Blunden era/. QOJ 3.7±0.4 
This Work 4.0 ±0.5 



The weak couplings for the electron are given by g e v = Tg — 
2Q e sin 9w(Q 2 ), andg^ = T e 3 . The axial integral receives sup- 
port from high Q 2 and we allow both a em and sin 2 9w to run. 
The dispersion relation that relates Im D^ z to Re D^ z is 



ReD 



rz 






dE' 



E' 



E' 2 - 



E 2—-yz(E')- 



rlmD^ 



(7) 



Blunden et al. II 1 Oil obtained axial results of the same order 
of magnitude as Re D^ z . Repeating a similar analysis we have 
also calculated ReD^ z . The two results for the axial contri- 
bution at the Q weil k energy are reported in Table |TT] As with 
the ReD!L calculation, differences between the axial results 
occur because of the structure function treatment. 

The goal of our paper is to describe our modifications to the 
electromagnetic structure functions. In Section II we present 
the steps taken to obtain Fj 2 m the resonance region. We fo- 
cus attention on this region since most of the support for the 
vector yZ-box integral comes from low Q 2 . These steps were 
not described in detail in our previous work and will allow a 
more thorough assessment of our Re D!L calculation. In Sec- 
tion III we describe an alternative modification for obtaining 
fJ 2 m the resonance region. This modification is similar to 
the one used by Gorchtein et al. [6jl and the close agreement 
to our original Re OY, z result suggests both modifications are 
equally valid, at least for the Q wea k kinematics. In Section 

IV we present our calculation of Fj in the resonance region 

yZ 
that parallels the analysis of Section II. We compare Fj and 

Re Oy Z values to those obtained by Blunden et al. Concluding 

remarks are contained in Section V. 



Following the normalization of the Particle Data Group |1J 
the resonant parts of these structure functions can be ex- 
pressed as a product of the polarization vector, £+ = 
l/\/2(0, — 1,— z',0), and hadronic tensors: 



7 yy(yz) 



JV-5-re.v 



ef«v (rZ) 
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1 J 



x(res,X\e+-J r (0)\N,s), 



(10) 



where X and s are the spin projections of the resonance and 
nucleon, respectively, and y (Z, V) is the electromagnetic (neu- 
tral vector) current. The factor of 2 is present in yZ-exchange 
to account for the different orderings. 

The above amplitudes can be evaluated by considering 
e + ■ J as a quark operator embedded between SU(6) wave 
function representations of the nucleon and resonances IU3I1 . 
This operator ignores the spatial wave functions, y/, and acts 
only on the flavor, 0, and spin, %, wave functions. Because the 
colorless portion of the total hadronic wave function is sym- 
metric, we are free to operate only on the third components of 
<p and x an d multiply the result by three. The amplitude can 
be expressed as 

(res,x\e + -j^\N 7 s) = 3xe[ 3 \ g f ) ) 

X (Yres(l> r esX\\u(k',X')e + ■ Yu{k,s')\WN<$>N%s), 

(11) 

where k {k') and s' {X 1 ) are the initial (final) momentum and 
spin projection for the struck quark. The superscript 3 over the 
quark electromagnetic and weak vector couplings, e q and g q v , 
indicates that the operators are acting only on the third quark. 
Using unit normalized quark spinors, 
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(12) 
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MODIFICATION OF THE STRUCTURE FUNCTIONS 

(x,Q 2 ) -> F^(x,Q 2 ) IN THE RESONANCE REGION. 



1,2V 



In our previous work we modified the Christy and Bosted 
fit Hill in the the resonance region Q 2 < 8 GeV 2 and W < 2.5 
GeV. Their fits for F^ r , o T , and O/, account for the contri- 
butions of seven resonances as well as a smooth background. 
Their description and computer code for their fit allowed us to 
separately modify the resonances and the background. For the 
resonances, our goal was to find prefactors for the summation 



Fj A =Y^C res xF } 



77 1 



(8) 



such that 



C re 



7 yZ 
1 
7 77 



(9) 



and choosing a frame where the gauge boson is propagating 
in the z-direction, the current reduces to 



u(k',X')e + ■yu(k,s') 



^ 



■q z S + ]&, (13) 



where m q is the constituent quark mass, P + = k\ + ikj, S + = 
l/2(oi + r(72), ^-is the momentum of the boson, and % are the 
usual two spinors. The Wigner-Eckart Theorem allows us to 
calculate a matrix element of P + as a constant times a matrix 
element of L + . 

After absorbing the spatial and momentum information, as 
well as the quark mass coefficient, into parameters A and B, 
Eq.dTOb becomes 
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Nitres 

(3) 



3x^(2^ (3) ) 



[Wn<PnXs\[AL+ +B5+] f I VrestyresXk) 



J?) 



X 3eq (YresQresXX \ [ AL + + B S+] \ WnQnXs) ■ 4 ) 



To calculate these amplitudes we operated the Hamiltonian on 
the SI/ (6) spatial (y/), flavor (<f>), and spin (%) wave functions 
of protons and resonances described by Close IU3I1 . As exam- 
ples, the proton and £>i3(1520) resonance are members of the 
(2 8 ,56) and (2 8 ,70) multiplets respectively and can be written 
as 



|2*,56> 



^*U-o(^ 



%S z =±l/2' 



.^A X M,A 



S z =±l/2, 

(15) 



|2 8 ,70)= £ <-/ = 3/2 7 z |LL z ,SSz) 

Jz=S z +L z 



5 [vffi (**•**£* -^ 



M.S 
z 



M,S V MA 
ls z 



r A x^ s 



(16) 



M,(A)S indicates a wave function with two elements that are 
(anti)symmetric. 

Inserting the Hamiltonian into the proton to Di3(1520) he- 
licity amplitudes gives 



A 1{L) - 3 x e (3) 


o# 3 


0(v«^^+i/ 
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[AL+ +BS+] \\I/ n <I>nXs 


= Te(- AlQ 


««(#)-«*(*£)" 




-V2B 10 


T j e u {gv) + -e d gv) 


) (17) 


and 


A I-3/2 = 3 X 4 3) (Sv 3) ) (Wres^resX+3/2 [AL+ + BS+] WnQnXs 


— ^ 

x/2 


10 


«> 


i{&)-eM\ 


3 


(18) 



where X is the spin projection of the resonance along the di- 
rection of the gauge boson momentum. The subscripts of A\o 
and Bio indicate the angular momentum dependence of the 
resonance's wave function. 

Obtaining Aio and Bio without relying on hadronic wave 
function requires additional phenomenological information. 
Data for both of theDi 3 (1520) and F J5 ( 1680) resonances (H 
U5I1 show that the polarization ratio 



A = 



I4/ 2 I 2 - 



1^3/2 1 



l4 /2 l 2 - 



■I4 /2 I 2 



(19) 



is close to — 1 for photoproduction, and changes to be close to 
+ 1 at higher Q 2 as the Aj,, amplitude dominates (in accord 
with perturbative QCD). Looking at the expressions for the 
Di3(1520), we conclude that 



A 10 (Q 2 = 0) = -V2B W (Q 2 =0) 



(20) 



\y 



and expecting A ' ,„ to dominate by a power of Q at high Q 



'1/2 

we choose a form with the correct limits 

Aio(e 2 ) 



Bw(Q 2 ) 



= -V2f l (Q 2 ) = -V2 



1 



1+G7A? 



(21) 



This leads to 



C Dl3 - 



(* 



■/i)(l-/i)- 



•3/? 



(Wi) 2 + 3/i 2 



n p.LO 



(22) 



where £>w = 1 — 4 sin 2 0w(O). A parallel analysis gives 



C Fl5 



1(1-/2) 
(l-/ 2 )2+2/ 2 2 



n p,LO 



(23) 



We used A 2 = A?; = 0.2 GeV 2 J5[]. These values were chosen 
by comparing our fits constructed using Close's analysis with 
MAID fits 11611 to transition amplitudes. Better agreement 
can be obtained by setting A 2 = 0.256 GeV 2 and A 2 = 0.635 
GeV 2 , but this more thorough analysis does not change the 
overall Re D^ z result by more than half a percent. 

Table [TIT] summarizes the prefactors for each resonance in 
the Christy and Bosted fit. The prefactor for the Roper reso- 
nance is easily calculated and is £) 2 -independent. The two Sn 
states belong to the same SU(6) multiplet as the Di3(1520), so 
Aio and Bio are the same for all three states, for valid SU(6) 
symmetry. The Si 1 states can mix. We have written above the 
results for the unmixed case. The unmixed Jp amplitude for 
the Si 1 (1650) is zero when the values of the quark charges are 
inserted; this is the Moorhouse selection rule ltl7H . If we ne- 
glect this amplitude also for the Z-boson case, the amplitude 
listed for the Sn (1535) gives a ratio 



c Sl 



-2/i 



l+2/i 



n p,w 



(24) 



Electroproduction of the Si 1 (1650) occurs because of mixing 
with the bare Sn(1535), and the above ratio is the same for 
both the Sn's. We have checked that including mixing makes 
little numerical difference. 

C res for 7 = 3/2 resonances are calculated by considering 
only the A/ = 1 portion of the current. This term is propor- 
tional to (e u — e c [). By substituting vector charges, C res for 
7 = 3/2 resonances is found to be (1 + Q 1 ^ ). 

The Christy-Bosted fit lies within 3% of nearly all data 
points. We estimated our modifications increased the uncer- 
tainty to 10%. 

The Christy-Bosted fit also accounts for a smooth back- 
ground. To model the yZ-box background we considered two 
limiting cases. In the low x limit, the light quark distributions 
are expected to be equal and the corrective coefficient is 



Cbkgd\. 



l^q=u,d,s ^ e qgyfq W 
2-iq=u,d,s{ e q) Jq( x ) 



1 



n p,LO 



(25) 



In the limit where there are only valence quarks 



r I _ 'Lq=ii,ii,d^ e q8vfq( x ) _ 2 n p,LO n ^ 

^■bkgd I valence quarks — ^ / \2 f ( \ -} x~-W ' ^ ' 

Lq=u,u,d\ e q) Jq\ x ) 3 

We used these limits as error bounds and their average as the 
background correction. Approximately half of the total con- 
tribution to Re Dy Z from the Christy-Bosted fit is due to this 
background modification. 
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FJ' is related to F' r by 



F yy - 



el 

p-q 



1 + — 



F™ 



M 2 e 2 



(27) 



yZ yZ 

We substituted F{ into the above expression to obtain F, . 

We also assumed the modifications were the same for both the 

transverse and longitudinal cross sections. 

Bosted and Christy iflSfl also have a fit for deuteron and 

neutron electromagnetic data which we used to modify the 

deuteron structure functions in [J9J] . The corrective ratios for 

the deuteron resonances are listed in Table [HI] Following the 

above analysis for the proton background, the limits to the 

deuteron background are 1 



e^ LO and4/5 



n p.LO 



0.0070 



0.0065 



0.0060 



tti 

^n 0.0055 



<x 



0.0050 



0.0045 



0.0040 
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III. 

x,Q 2 ) 



ALTERNATIVE MODIFICATION OF 

IN THE RESONANCE REGION 



■Ff 2 ( X) Q^ 



1,2 V 



The corrective prefactors for the Christy-Bosted fit can be 
modeled using a different technique. The vector contribution 
to the Z-boson transition amplitudes can be isospin rotated 
into a sum of electromagnetic transition amplitudes, py —t N* 
and ny — > N*. Neglecting strange quark contributions, these 
amplitudes are 



FIG. 1 : Re D^ z as a function of incoming electron energy. The black 
curve is the result from our previous work and uses helicity ampli- 
tudes given by the quark model. The blue, dot dashed curve is the 
result using corrective ratios from the PDG The red, dashed line is 
the result from using corrective ratios constructed with MAID he- 
licity amplitudes. All three models use the same modifications for 
isospin 3/2 resonances and the smooth background. 



T r(zy) 



(N;\J^ v >\p) = e u (g u v)(N;\uYnu\p) + e d ( g ^)(N;\dy fl d\p) 



and 



(28) 



(N:\4\n) = e u {K\UYiiu\n) +e d (N:\dy fi d\n). (29) 



After performing an isopin rotation the neutron amplitude be- 
comes 



{N* n \jJ l \n}=e ll (N* p \dy fl d\p) + e d (N* p \uy^u\p) 
Further algebra on these amplitudes reveals 



(30) 



(K\j z /\p) 



i 



(1 -4sin 2 0w(O))(Ag^|p} - AK\Jl\n) 



(31) 



C re ? can now be written as 



^res — -^ 



U(aV) 2 



n P,LO _ 



i y,p a r-« 



U(aV) 2 



(32) 



Gorchtein et al. la] constructed their C res expressions us- 
ing photoproduction amplitudes listed in the Particle Data 
Group [12]. Thus, their corrective prefactors lack Q 2 - 
dependence. To account for the amplitudes' Q 2 -dependence, 
fits from MAID iflal can also be used. 

Fig.[T]shows Re D!L calculated using both the quark model 
and MAID treatments of the structure functions. Better agree- 
ment between MAID and the quark model was naively ex- 
pected as the MAID fits were used to parameterize A 2 2 . The 



overall smaller value for Re D^ z calculated by MAID is al- 
most entirely due to the Roper resonance. For the Roper, the 
quark model calculates a constant corrective prefactors while 
the MAID ratio rapidly approaches Q^ as Q 2 increases. 
The differences in the Roper resonance corrective prefactors 
were also the primarycause for the different deuteron asym- 
metry predictions in H . 

Another notable feature of Fig. Q] is that Re DY, Z hardly 
changes when the corrective ratios are calculated using PDG 
photoproduction amplitudes in place of the Q 2 -dependent 
quark model. ReD*^ calculated using the quark model also 
remains relatively unchanged when using different values for 
Aj 2 values. Both features are due to low Q 2 values domi- 
nating the integral. Indeed, an analysis of the integral indi- 
cates that the mean Q 2 value is 0.4 GeV 2 . In applications with 
higher Q 2 , such as the calculation of the deuteron asymmetry 
in Q9Q, the quark model and photoproduction corrective pref- 
actors give quite different values. 



It is important to note that Gorchtein et al. la] do not use 
the Christy-Bosted background in their analysis. For the back- 
ground they instead use the average of two Generalized Vector 
Dominance (GVD) models ll9l 12011 . isospin rotated for appli- 
cation to the yZ-box and extrapolated down to the resonance 
region. This averaging is the largest source of uncertainty for 
the Gorchtein et al. calculation. Recently, it has been claimed 
that this background uncertainty has been overestimated iQ]. 



TABLE III: The seven Christy-Bosted resonances along with their electromagnetic helicity amplitudes along and corresponding corrective 
prefactors for both the proton and deuteron. The (pZ —¥ Nl) helicity amplitudes are calculated by substituting e q — > g v = T 3 — 2e ? sin 9w . 
The (ny— > N*) and (nZ —t N*) helicity amplitudes are calculated by exchanging e„ ++ e d and g v ++ g^, respectively, in the proton analysis. 
The corrective prefactor for the background is also included. 
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proton electroproduction amplitudes 
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A 3/2 = ^ A 2o( 2e « + ^) 
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IV. MODIFICATION OF STRUCTURE FUNCTION 

x, Q 2 ) -+ F^(x, Q 2 ) AND THE CALCULATION OF 



F rr 



ReD^ 



Blunden et al. 11011 split their Re Dz z analysis into elastic 

(W 2 = M 2 ), resonance (W 2 < W 2 < 4 GeV 2 ), and deep in- 
elastic scaling (W 2 > 4 GeV 2 ) regions. To allow for an easier 
comparison between our analysis and theirs, we used the same 
energy regions. 

As previously mentioned, the average Q 2 value within the 
Re DY, Z integral is about 0.4 GeV 2 . In contrast, the average Q 2 

value within the ReQL integral is about 80 GeV 2 . Thus, the 
axial contribution to the yZ-box diagram is less sensitive to 
the modifications of the structure functions in the resonance 
region. Because the axial box integral, Eq.((5), receives strong 
support from high Q 2 , we follow the example of Blunden et 
al. and use one loop running values of a (Q 2 ) and sin" 0w {Q 2 ) 
in its evaluation. 

In the scaling region Fj can be directly calculated using 
parton distribution functions 

?7Z 



For Q 2 < 1 GeV 2 and W 2 > 4, we used Model I modifica- 
tion to the PDFs discussed by Blunden et al., with A 2 = 0.7 
GeV 2 and Qq = 1 GeV 2 . We approximated this uncertainty 
to be 10%. For the elastic contribution, we also follow the 
technique used by Blunden et al. lUOll . 

The most significant departure from the Blunden et al. anal- 
ysis is in the resonance region. In this region Blunden et al. 
constructed Fj using axial current parameters of Lalakulich 
et al. 12311 . Lalakulich et al. obtained their parameters through 
a PCAC analysis of pionic decays of baryons. Their fit ac- 
counts for four resonances but makes no attempt at estimating 
a smooth background, defering the determination of its form 
to future experiments. As an aside, Lattice QCD calculations 
have reached a sufficient level of accuracy to calculate axial 
form factors JUIH. 

Instead of repeating the Blunden et al. resonance region 
analysis, we constructed Fj by once again modifying the 
Christy-Bosted fit. Not only does this modification provide 
a smooth background, it also accounts for three more reso- 
nances. In our analysis of the resonance region we repeated 
the technique outlined in Sec. [II] In the non-relativistic limit, 



Ff(x, Q 2 ) = Y,2e q g 9 A (q(x, Q 2 ) - q(x, Q 2 )) . (33) 1*1 « m q , the axial current becomes 



Blunden et al. use PDFs from 12 1H . We chose PDFs given 
by CTEQ [22]. To be conservative we considered a 10% 
uncertainty for this fit. 



ii(k',X')e + ■ yy 5 u{k,s') = V2tf,S + i,,. (34) 



Continuing the use of the parameters in Sec. [II] Fj can be 
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FIG. 2: Elastic (black, solid curve), resonance (red, dashed curve), 
and model I (blue, dot dashed curve) contributions to the axial box. 
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FIG. 3: The axial box. We also add the axial box to our previous 
vector calculation [50 to obtain the total box. 



expressed as 



V. CONCLUSIONS 
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WresQresXl) 



X 3e^ ' (yfresQresXX \ [ AL + + B S+] \ WPnXs) , (35) 

where v is the energy of the exchanged boson. For our cal- 
culation we took the mass of the struck quark m q to be 0.3 

GeV. Table II V I summarizes the corrective prefactors to obtain 

yZ. yZ 

F 3 ' . As with the corrective prefactors for F i ',, we estimate 

the uncertainty of the Fj prefactors to be 10%. 

The smooth background is once again modified by taking 
the low x and valence quark limits. For low x, a quark and 
anti-quark are equally likely to be struck. Thus, 



Cbkgd |jc-M) — 



I.q=u ,d,.s 2e qglfq( x ) 

u 1 d,s\ e q)~"Jq\ x ) 



2H= 



= 0. 



(36) 



In the limit where valence quarks are equally likely to be 
struck 



^bkgci I valence quarks 



Hq=u,uJ^qg q A f q (x) _ 10 
2^q=u 1 u,d\ e q) Jq\ x ) * 



These limits were taken as the uncertainty bounds and their 
average as the modification for the smooth background. 

We also calculated f7 for the deuteron in Q. The correc- 
tive ratios for the deuteron resonances are listed in Table |IV] 
Following the above analysis for the proton background, the 
limits to the deuteron background are and 18/5. 

Figs.|2]and [3]display the results for Re Qj z . As can be seen, 

the scaling region dominates. At the Q wea k energy, ReQL = 
0.0040 ±0.0005. 



Adding the axial box to our original vector box calcula- 
tion J5D, our constituent quark model yields a total yZ-box 
value of 



ReD rZ (£ = 1.165 GeV)| to tai = (9.7 ± 1.4) x 10" 



(38) 



The errors from both the axial and vector calculations were 
added directly. If added in quadrature, the uncertainty reduces 
to 1 x 10~ 3 . 

The total yZ-box value from Blunden et al. lUOIl is 



ReD rZ (£ = 1.165 GeV)| total = (8.4+^) x lO" 3 . 



(39) 



These two calculations are in agreement within uncertain- 
ties. Each calculation also has error bounds below the error 
budget of the Q wec ,k experiment. 

The question remains which calculations the Q wea k collabo- 
ration should use in their analysis. The disagreement between 
the various calculations is largely due to the treatment of the 
yZ structure functions in the resonance region. We believe 
the collaboration will be equally well-served by either Re D^ z 
calculation. Re D^ z is not very sensitive to the resonance re- 
gion modifications since its integrals get much of their support 
(37) from high Q 2 . Fj in the scaling region can be constructed us- 
ing fits to parton distribution data. 

Which ReD^ z calculation to use is more open to debate. 
The vector integrals receive much of their support from the 
resonance region and are thus sensitive to the modification 



^12 "~ * F\2 ■ I n S ec - Ell we showed that there is little differ- 
ence between modifying the Christy-Bosted resonance fits us- 
ing our constituent quark model [ 5 ] or photoproduction am- 
plitudes from the Particle Data Group (as in [j6j]). Differ- 
ences arise between |5|] and la] because of the treatments of 
the resonance region background. We continue modifying the 
Christy-Bosted background fit while Gorchtein et al. modify 



TABLE IV: The seven Christy-Bosted resonances along with their axial helicity amplitudes and corrective prefactors for both the proton and 
deuteron. The neutron amplitude is calculated by exchanging g u A <-> g^. 
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two GVD fits to low Q 2 , high W 2 data and extrapolate them 
down to the resonance region. We believe our modification is 
more satsifactory since it does not involve any extrapolations. 
We cannot comment on the vector calculation of |Ql since they 
provide few details of their model. 
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